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The current -voltage {I-V) characteristics of single-electron transistors (SETs) have been mea- 
sured in various electromagnetic environments. Some SETs were biased with one-dimensional ar- 
rays of dc superconducting quantum interference devices (SQUIDs). The purpose was to provide the 
SETs with a magnetic-field-tunable environment in the superconducting state, and a high-impedance 
environment in the normal state. The comparison of SETs with SQUID arrays and those without 
arrays in the normal state confirmed that the effective charging energy of SETs in the normal state 
becomes larger in the high-impedance environment, as expected theoretically. In SETs with SQUID 
arrays in the superconducting state, as the zero-bias resistance of the SQUID arrays was increased 
to be much larger than the quantum resistance Rk = h/e 2 « 26 kf2, a sharp Coulomb blockade 
was induced, and the current modulation by the gate-induced charge was changed from e periodic 
to 2e periodic at a bias point < |V| < 2Ao/e, where Ao is the superconducting energy gap. The 
author discusses the Coulomb blockade and its dependence on the gate-induced charge in terms of 
the single Josephson junction with gate-tunable junction capacitance. 
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I. INTRODUCTION 

The single-electron transistor (SET), which consists of 
two small-capacitance tunnel junctions in series and a 
gate electrode, is an important example of the single- 
electron-tunneling circuit. Although the SET is more 
complex than the single junction, the SET has an advan- 
tage over the single junction in that the single-electron 
charging effects typified by "Coulomb blockade" are eas- 
ily observed in the SET. As a result, the SET can be 
used, e.g., as an electrometer with an extremely high sen- 
sitivity, <§; e/VHz (e.g., Ref.QJ). When the electrodes of 
the tunnel junctions in the SET are in the superconduct- 
ing state, the SET is also important from the viewpoint 
of quantum computing; the first one-quantum-bit opera- 
tion in solid-state electronic devices was demonstrated in 
a circuit with small-capacitance Josephson tunnel junc- 
tions, and the state of the quantum bit was read out 
by measuring the Josephson-quasiparticle (JQP) current 
of a superconducting SET^ The supercurrent of super- 
conducting SETs has also been studied extensively. 3 One 
expects that the supercurrent is 2e periodic in the gate- 
induced charge. In many experiments, however, a period 
of e has been seen, which suggests the existence of sub- 
gap quasiparticle states^ or "quasiparticle poisoning." 
The periodicity of the current was also discussed at fi- 
nite voltages4^& The periodicity depends strongly on 
how the leads connected to the sample are filtered^ so 
that the "quasiparticle poisoning" would be viewed as an 
environmental effect. 

Theoretically^ the details of single-electron charging 
effects in SETs depend on the impedance of the electro- 
magnetic environment, Z(u>), both in the normal state 
and in the superconducting state. It is predicted that as 



Re[Z(w)] is increased, the voltage scale of the Coulomb 
blockade in the normal-state current-voltage (I-V) char- 
acteristics becomes larger, and in the superconducting 
state, the supercurrent is replaced by the Coulomb block- 
ade. The experiments^^ on superconducting SETs 
mentioned above were done in the low-impedance en- 
vironment, where Re[Z(o;)] was much smaller than the 
quantum resistances: Rq = h/(2e) 2 ss 6.5 kfl for Cooper 
pairs and Rk = h/e 2 ~ 26 kfi for quasiparticles. In or- 
der to obtain higher Re[Z(w)], thin-film resistors^ can 
be used for the on-chip leads. A superconducting SET 
was biased with «50 ktt Cr resistors, and a Coulomb 
blockade was observed in the I-V curved Similar re- 
sistors (2 — 20 kfi) were also employed for the study of 
Cooper-pair cotunneling in superconducting SETs^i A 
drawback of these thin-film resistors is that the resistance 
is not tunable, i.e., the same SET cannot be measured 
in different environments. A tunable environment has 
been realized by capacitively coupling a two-dimensional 
electron gas to a SET. 12,13 In the work presented here, 
a simpler method is employed to create tunable environ- 
ments for superconducting SETs. 

We use one-dimensional (ID) arrays of dc supercon- 
ducting quantum interference devices (SQUIDs) for the 
on-chip leadsji^i^ which can be fabricated simultane- 
ously with superconducting SETs. The SQUID configu- 
ration enables us to vary in situ the effective impedance 
of the arrays by applying a weak external magnetic field 
(1 — 10 mT) perpendicular to the SQUID loop. The su- 
perconducting SET in our samples, on the other hand, 
does not have a SQUID configuration, and therefore its 
parameters are practically independent of the external 
magnetic field. The zero-bias resistance of SQUID arrays 
at low temperatures can be controlled over several orders 



of magnitude (e.g., Fig. 13 of Ref . Il5|) . and SQUID arrays 
are especially suitable for achieving a high-impedance 
environment. In fact, distinct Coulomb blockade was 
observed in single Josephson junctions by biasing with 
SQUID arrays^ When the sample is driven to the nor- 
mal state, the array leads are no longer tunable, but the 
resistance can be much larger than Rk- Hence, it is 
still possible to study the environmental effects on SETs 
by comparing the samples with and without array leads. 
In this paper, the environmental effects on SETs is dis- 
cussed both in the normal state and in the superconduct- 
ing state. An emphasis is placed on the case of a high- 
impedance environment, which has not been thoroughly 
experimentally studied. 

II. EXPERIMENT 

A. Fabrication of small-capacitance tunnel 
junctions 

The tunnel junctions (AI/AI2O3/AI) were fabricated 
on a SiC-2/Si substrate with Au/Ni bonding pads. We 
employed a process based on electron-beam lithography 
and double-angle shadow evaporation, which is similar to 
the one described in Ref. The evaporation of Al was 
done at an rate of 0.1 — 0.2 nm/s in a vacuum system with 
the base pressure of < 10 -8 Pa. During the evaporation, 
the pressure was usually (2 — 4) x 10~ 6 Pa. The AI2O3 
tunnel barrier was formed by exposing the base Al layer 
to 1 — 20 Pa of O2 for 0.5 — 2 min. before the deposition of 
the top Al layer. The thickness of the barrier determines 
the normal-state resistance of the junction per unit junc- 
tion area. The normal-state resistance is a key parameter 
of the sample, and will be discussed in Sec. Ill CI 

Figure ^ shows a scanning electron micrograph of a 
sample. The SET is biased with two pairs of leads en- 
abling four-point measurements. A part of each lead close 
to the SET consists of a ID array of dc SQUIDs. 

B. Measurement at low temperatures 

The samples were measured in a 3 He- He dilution re- 
frigerator (Oxford Instruments, Kelvinox 400) mainly at 
T = 0.02 — 0.6 K, and the normal-state resistance was 
determined at T = 1.8-4.5 K. (Note that the su- 
perconducting transition temperature of Al is 1.2 K.) 
The temperature was calculated from the resistance 
of a ruthenium-oxide thermometer fixed to the mixing 
chamber. Magnetic fields on the order of 1 — 10 mT 
were applied by means of a superconducting solenoid. 
In this magnetic-field range, the temperature error of 
ruthenium-oxide thermometers due to magnetoresistance 
is negligibly small (e.g., a typical value of the error at 
T = 0.05 K is less than 0.1%) The samples were placed 
inside a copper can which was thermally connected to the 
mixing chamber. Because there was no low-temperature 
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FIG. 1: (Color online) A scanning electron micrograph of an 
AI/AI2O3/AI single-electron transistor (SET) with SQUID- 
array leads. An equivalent circuit for the SET is shown in the 
top left corner. 



filtering, we inserted low-pass filters between the cables 
connected to the cryostat and the measurement circuit. 

The I—V curve of the SET was measured in a four- 
point configuration (see Fig. QJ. The bias was applied 
through one-pair of leads, and the potential difference 
was measured through the other pair of leads with a 
differential voltage amplifier (DL Instruments, 1201, > 
1 Gfi input impedance). The current was measured with 
a current amplifier (DL Instruments, 1211). The SQUID 
arrays could be measured in a two-point configuration 
(same current and voltage leads). When we measure 
the arrays on the same side of the SET (e.g., 1+ and 
V+ in Fig. nj, current does not flow through the SET, 
and the series resistance of two arrays is obtained. The 
zero-bias resistance measured in this configuration will 
be discussed in Sec. Ill CI 



C. Characterization of the samples 

In this work, we measured the three pairs of SETs 
listed in Table |U Each pair (e.g., samples la and lb) 
was fabricated simultaneously on the same chip. One 
SET has SQUID-array leads with N = 65 junction pairs 
in each array, while the other SET has no SQUID-array 
leads (N = 0). For all the samples, the junction area was 
designed to be 0.1x0.1 ^m 2 in the SET and 0.3x0.1 /mi 2 
in the SQUID arrays. 

The uniformity of the SQUID arrays could be esti- 
mated by measuring the normal-state resistance in all 
two-point configurations (six in total). In some config- 
urations, current flows through the SET and thus the 
normal-state resistance R n of the SET has to be sub- 
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TABLE I: List of the samples. R„ is the normal-state resis- 
tance of the single-electron transistor, r' n is the normal-state 
resistance per junction pair of the SQUID-array leads, and C g 
is the capacitance between the island electrode and the gate 
electrode. Samples lb, 2b, and 3b do not have SQUID arrays 
in the leads. 



pendicular to the SQUID loop, 



Sample 


Rn (kfi) 


r' n (kfi) 


C g (aF) 


la 


82 




6.0 


lb 


102 


5.6 


6.0 


2a 


57 




4.8 


2b 


85 


4.1 


4.9 


3a 


17 




4.8 


3b 


23 


1.4 


4.8 
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FIG. 2: (Color online) Zero-bias resistance of two SQUID- 
array leads connected in series on the same side of the SET 
vs external magnetic field for sample 2b at T = 0.02, 0.3, and 
0.6 K. 



tracted from the measured resistance in order to obtain 
the series resistance of two arrays. For our samples, the 
results of the two-point measurements agreed with one 
another within 3%-4%. By averaging all the results, the 
normal-state resistance r' n per junction pair of the arrays 
is calculated, and shown in Table [I] If the tunnel-barrier 
thickness is identical for all junctions on the same chip, 
one would expect R n /r' n w 12 from the junction area. 
From Table Q] the ratio is 16 — 21, which is the correct 
order of magnitude. 

From the normal-state resistance, I calculate the 
Josephson energy, which is an important parameter of 
the samples. For example, the maximum Josephson en- 
ergy E' J0 between adjacent islands in the SQUID arrays 
is given by the Ambegaokar-Baratoff formula, 



E' 



.70 



8e 2 r' 



(1) 



where Ao is the superconducting energy gap (« 0.2 meV 
for Al). Because of the SQUID configuration, the ef- 
fective Josephson energy E'j in the arrays is modulated 
periodically by applying an external magnetic field B per- 



E'j = E' 



JO 



BA 



(2) 



so long as B is sufficiently smaller than the critical field, 
where $0 = h/2e — 2 x 1CU 15 Wb is the superconducting 
flux quantum, and A is the effective area of the SQUID 
loop, which is 0.7x0.2 /Ltm 2 in our samples. Equation @ 
is a key to understand Fig. [2 which shows that the zero- 
bias resistance R' of two SQUID arrays in series oscil- 
lates as a function of B until B becomes comparable to 
the critical field. Note that the first peak appears at 
B = 7 mT, where the normalized flux BAJJ&Q is 0.5 for 
our samples with A = 0.14 ^m 2 . In Ref.El the SQUID 
array was modeled as a network of capacitors and in- 
ductors, and the real part of the array impedance was 
shown to be proportional to 1/ y/E'j at low enough fre- 
quencies. At B > 60 mT in Fig. [21 the superconductivity 
is suppressed, and R' takes a value comparable to 2Nr' n . 
The small temperature dependence in this magnetic-field 
regime is due to the single-electron charging effect. In 
general, the I-V curve of SQUID arrays is nonlinear (e.g., 
Fig. 3 of Ref.El, and therefore, SQUID arrays are not 
described by a linear impedance models However, it 
would be a good approximation^ to use R' Q for charac- 
terizing the electromagnetic environment of the SET. I 
also assume that the SET is in a high-impedance envi- 
ronment when R' 3> Rk- As for the samples without 
SQUID arrays, the resistance of the leads is always much 
smaller than Rk, and the SET is in a low-impedance 
environment. The Josephson energy in the SET is calcu- 
lated in a similar way, but note that it does not oscillate 
as a function of B because the junctions in the SET do 
not have a SQUID configuration. 

Another important parameter is the charging energy, 
which is inversely proportional to the capacitance of the 
junction. The capacitance can be estimated from the 
junction area with a specific capacitance on the order 
of 10 2 fF//^m 2 . In our samples, the junction area is 
much larger, or the charging effect is much weaker, in 
the SQUID arrays than in the SET. In the normal state, 
for example, the I-V curve of the SQUID arrays is al- 
most linear (data not shown, similar to the bottom curve 
111 Pig.3ofRef.Ii3) even when the SET shows a well de- 
veloped Coulomb blockade like in Fig.HJb), which would 
be favorable when the arrays are used as leads. 

The capacitance C g between the island electrode and 
the gate electrode is determined by the period of the gate 
modulation in the normal state. An example is shown in 
Fig. 13 Though the design of the SET is the same for all 
the samples, the values of C g in Tableware divided into 
two groups: 6 aF for samples la and lb, and <5 aF for 
the others. The reason could be that we used two SiC^/Si 
wafers whose thickness t of the SiC>2 layer is different, and 
that we fabricated samples la and lb on one wafer (t = 
0.2 fxm) and samples 2a-3b on the other (t = 0.5 fim). 
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FIG. 3: (Color online) Modulation of the current by the gate 
voltage in the normal state (B = 0.1 T) at T = 0.02 K for 
sample la. The bias voltage is 0.04 mV. 
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FIG. 4: (Color online) Double-junction system connected to 
a voltage source V a via electromagnetic environment Z{ui). 



III. RESULTS AND DISCUSSION 

A. Single-electron transistors in the normal state 

At T = 0, the I-V curve of the SET in the normal 
state is expected to have a region where I — 0. Such a 
region is called a "Coulomb gap," and depends not only 
on the gate voltage V g but also on the electromagnetic 
environment £ The circuit considered in the theory^ is 
shown in Fig.^J where Ci and Ri (i — 1,2) are the capac- 
itance and the tunnel resistance, respectively, for the ith 
junction, and Z(lo) is the impedance of the electromag- 
netic environment. The theory assumes that Ri 3> Rk 
and Ri ^S> Re[Z(ui)]. In Fig. the gate is not drawn 
explicitly because when d 3> C g , which is the case in 
our samples, the influence of the gate can be included in 
the effective island charge, 



q = ne + C„V g + Q , 



(3) 



where n is an integer and Qq is the background charge. 

In the theory, the drift of Qo is not considered, and 
thus for simplicity, Qq = in Eq. @. In the experiment, 
however, the drift often becomes a problem. In sam- 
ples la, lb, 3a, and 3b, the drift was usually < O.le/day, 
which is negligible for the discussion in this paper. In 
samples 2a and 2b, on the other hand, the drift was not 
slow enough (~ e/hour) compared to the time needed to 
measure hundreds of I-V curves at different values of V a 
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FIG. 5: The theoretical region of zero current at T — 
(Coulomb gap) for single-electron transistors with C2/C1 = 
1.5 in (a) the low-impedance environment and in (b) the high- 
impedance environment. The source voltage V s is in units of 
e/(Ci + Ci) and the gate voltage V g is in units of e/C g . 



in our experimental setup. Thus, for samples 2a and 2b, 
we mainly measured I vs V g keeping the source voltage 
V s constant. 

Figure El shows the theoretical Coulomb gap on the 
V g —V s plane for Z(u>) = (low-impedance environment) 
and for Re[Z(o;)] 3> Rk (high-impedance environment). 
The Coulomb gap is periodic in V g for both the cases. 
However, the magnitude of the gap is different, and e.g., 
at V g = 0, ±e/C g , ±2e/C g , the gap in the high- 
impedance environment is about twice as large as that in 
the low-impedance environment. It is important to note 
a qualitative difference at V g = ±0.5e/C g , ±1.5e/C g , 
±2.5e/C ff , .... The gap vanishes in the low-impedance 
environment, but survives in the high- impedance envi- 
ronment. These theoretical predictions have been con- 
firmed in our experiments. Figure HO shows the normal- 
state I-V curves at T = 0.02 K for samples la and lb. 
We drove the samples into the normal state by applying 
a magnetic field of 0.1 T perpendicular to the substrate. 
In the normal state, the SQUID arrays of sample lb have 
= 1.4 MO > R K at T = 0.02 K, and thus, I assume 
that the SET in sample lb is in a high-impedance envi- 
ronment. The SET in sample la, on the other hand, is 
in a low-impedance environment because it is not biased 
with the arrays. In the lower data set (q = 0.0), the size 
of the Coulomb gap in Fig. HJb) is about twice as large 
as that in Fig. Eta). Moreover, in the upper data set 
(q = 0.5), the I-V curve in Fig. Et a ) is linear, i.e., no 
Coulomb gap, while the curve in Fig. Etb) still shows a 
considerable nonlinearity. 

I have also calculated the I-V curves based on the 
theory assuming R± + R2 = R n - For the lower data 
set in Fig. Etb), a perfect agreement is obtained with 
R1/R2 = 1.0, C1/C2 = 1.0, Ci + C 2 = 1 fF, and 
k B T/E c = 0.1 (T = 0.09 K), where Ec = e 2 /2(C 1 + C 2 ). 
These parameters are employed for all the solid curves in 
Fig. El The solid curve does not always reproduce the 
experimental data very well, but I nevertheless conclude 
that our experiments have demonstrated the environmen- 
tal effects predicted by the theory. I emphasize that the 
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FIG. 6: (Color online) Current-voltage characteristics of 
single-electron transistors in the normal state in (a) low- 
impedance environment for sample la and in (b) high- 
impedance environment for sample lb. The current and the 
voltage are in units of e/(Ri +7?2)(Ci + C2) and e/(Ci + C2), 
respectively. The solid circles are experimental data at B = 
0.1 T and T = 0.02 K, where R x +R 2 = Rn and C1+C2 = 1 fF 
are employed. The solid (dashed) curves are the theoreti- 
cal prediction for k B T / E c = 0.1 (0.3), R l /R 2 = 1.0, and 
C1/C2 = 1.0. The upper (lower) data sets are for q = 0.5 
(0.0). The origin of the current axis is offset for each curve 
for clarity. 
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FIG. 7: (Color online) Offset voltage defined in Eq. @ as 
an estimate of the Coulomb gap for (a) sample la and (b) 
sample lb. The broken lines are the theoretical predictions 
for (a) low-impedance environment and (b) high-impedance 
environment. The units of the axes are the same as in Fig. O 
and the parameters, Ri, R2, Ci, and C2, are the same as in 
Fig.© 

solid curves are calculated with the same set of parame- 
ters, and the discrepancy is reduced if the parameters are 
adjusted for the calculation for each data set. The dashed 
curves for the lower data set of Fig. Efa) and for the up- 
per data set of Fig. Efb), which agree better with the 
experimental data, are obtained by raising the tempera- 
ture in the calculation to k B T/E c = 0.3 (T = 0.28 K), 
keeping the other parameters constant. 

At finite temperatures, nonzero current flows even 
within the Coulomb gap. Therefore, it is not straightfor- 
ward to determine the Coulomb gap from the experiment. 
As an estimate of the Coulomb gap, I have calculated the 



offset voltage V Q g by fitting 

I<x(V-V oS ) (4) 

to the data in the high-bias regime, / > 0.3e/(i?i + 
i? 2 )(Ci+C 2 ) iorVaff > Oand/ < -0.3e/ (R 1 + R 2 )(C 1 + 
C2) for V s < 0. The results for samples la and lb are 
plotted in Fig.|7|together with the theoretical predictions 
for the low-impedance environment and high-impedance 
environment, respectively. Here, I used the same param- 
eters as in Fig. |SJ The experimental V ff is consistent 
with the theoretical prediction. 

In fact, one of the assumptions in the theory for the 
high-impedance case, Ri ^> Re[Z(oj)], is not fulfilled in 
the experiment. In the normal state, R n < R' in sam- 
ple lb, so that the SET in sample lb is current-biased. 
However, I am still convinced that the comparison with 
the theory is meaningful, because we measured the SET 
in four-terminal configuration, and because the measured 
I-V curves are qualitatively explained within the theory. 
Furthermore, the parameters used in the calculations for 
Fig.|H|are reasonable. It is common to observe that the ef- 
fective electron temperature becomes considerably higher 
than the temperature of the mixing chamber in the exper- 
iments of small-capacitance tunnel junctions^ The tem- 
perature difference is likely to be large when the cryostat 
leads are not filtered at low temperatures, which is the 
case in our cryostat. Most importantly, C\ + C2 = 1 fF is 
consistent with our junction size, and as I mentioned ear- 
lier, much larger than C g . This value of the total capac- 
itance will be used again in the discussion in Sec. IIII Bl 



B. Single-electron transistors in the 
superconducting state 

For the samples with SQUID-array leads, the same su- 
perconducting SET can be studied in different electro- 
magnetic environments by changing the external mag- 
netic field on the order of a few mT. In this work, I 
look at the low-bias region of \V\ < 2Aq, where the I-V 
curve is sensitive to the state of the electromagnetic en- 
vironment. In all the samples with SQUID-array leads, 
the I-V curve of the superconducting SET developed 
a Coulomb blockade as R' was increased by tuning the 
field. Figure [S] shows the I-V curves of the SET in sam- 
ple lb in two different environments. As I mentioned ear- 
lier, the parameters of the SET should be independent of 
the field, because the SET does not have a SQUID con- 
figuration and the field applied in Fig. [S] is much smaller 
than the critical field. The electromagnetic environment 
for the SET (the SQUID arrays), on the other hand, is 
strongly varied with the field as we have seen in Fig. [21 
The behavior of the SET demonstrated in Fig.[S]does not 
result from the magnetic-field influence on the I-V curve 
of the SET, but rather from an environmental effect on 
the SET. I note here that the SETs without SQUID ar- 
rays were also measured at B = and at B « 7 mT, and 
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FIG. 8: (Color online) Sets of the current-voltage curves 
of the same superconducting SET (sample lb) in different 
environments: (a) R' = 0.2 Mft (B = 0), and (b) R' = 
0.3 Gtt (B — 6.8 mT), where -Ro is the zero-bias resistance 
of two SQUID-array leads connected in series. Both in (a) 
and (b), from bottom to top, the normalized gate-induced 
charge (C g V g + Qo)/e increases from —2.0 to +2.0 in steps of 
0.5. The origin of the current axis is offset for each curve for 
clarity. All the curves were measured at T = 0.02 K. 



that the I-V curves were almost the same. The well- 
developed Coulomb blockade in Fig. 02b) indicates that 
the SET is in a high- impedance environment, which is 
consistent with R' Q = 0.3 G£l ^> Rq. From the viewpoint 
of R' , the SET is not in a low- impedance environment 
in Fig. |S[a) because R' Q = 0.2 Mfi is already larger than 
Rq. Thus, in order to study the low- impedance case, 
we have to measure samples without SQUID-array leads. 
Rather surprisingly, e.g., in sample la, the I-V curve and 
its dependence on the gate-induced charge were qualita- 
tively almost the same as Fig. Ufa) > an< ^ the supercur- 
rent was too small to be seen in the I-V curve. This is 
due to "quasiparticle poisoning," which is likely to occur 
in our cryostat which has no low-temperature noise fil- 
ters. In addition, our biasing scheme is not ideal for the 
supercurrent measurement because the samples without 
SQUID-array are voltage biased. The supercurrent was 
not detected in samples 2a or 3a, either. 

"Quasiparticle poisoning" is also suggested by the 
e periodicity in the current modulation by the gate- 
induced charge, which was observed in all samples with- 
out SQUID-array leads. The e periodicity was also found 
in sample lb when R' = 0.2 Mfi, as shown in Fig.|j5Ja), 




FIG. 9: (Color online) Modulation of the current by the gate- 
induced charge in the same SET as in Fig. [g] The potential 
drop across the SET is 0.15 mV both in (a) and (b). The 
state of the environment for the SET in (a) and (b) is the 
same as in Fig. [Hfa) and Fig. |Ub) , respectively. 



where the potential drop across the SET is 0.15 mV. 
When R' is increased to 0.3 Gfi [Fig. E^b)], it was re- 
placed by 2e periodicity, which indicates Cooper-pair 
transport. The SQUID arrays with sufficiently large R' 
acted as a filter that suppresses the contribution of quasi- 
particles to the charge transport. This change in the peri- 
odicity was also found in sample 2b. In sample 3b, how- 
ever, only e periodicity was seen probably because the 
maximum of R' was only 9 Mil. The curve in Fig. 0b) 
is calculated from the I-V curves, because we could not 
fix the potential drop across the SET when it is much 
smaller than that across the SQUID arrays in our exper- 
imental setup. 

Below I focus on the case of R' ^> Rq, and dis- 
cuss the Coulomb blockade of Cooper-pair tunneling. 
In the high- impedance environment, a current-biased 
single Josephson junction is expected theoretically^^ 
to have a "back-bending" I-V curve, which has been 
experimentally 814,15 confirmed. The "back-bending" is 
also seen in Fig. H^b) for certain values of the gate- 
induced charge. When the "back-bending" is clearly seen 
in the SET I-V curve, I define the blockade voltage Vb 
as the local voltage maximum (minimum) for Vb > 
(Vb < 0) in the low-current part of the I-V curve. The 
measured Vb is plotted in Fig. ^| as a function of the 
gate-induced charge. The dotted and broken curves rep- 
resent V at constant / calculated from the I-V curves. 
These curves would also characterize the Coulomb block- 
ade. The analytic expression of Vb at low temperatures 
has been obtained for single Josephson junctions as 



V b « 0.25 



C 



for Ej/E c < 1, and 



V b cx 



— ( —L\ 



3/4 



C \E C 



exp 



EcJ 



1/2 



(5) 



(6) 



for Ej j Eq 3> 1 , where C is the capacitance of the single 
junction. 



7 



0.2 

> 

a 


- 


■ i I i i i i I i i 


/\ 


^ 0.0 

u 






- 








■ - . - / 


-0.2 


V 


V 

. . 1 .... 1 . . 






-2 -1 


1 


2 



FIG. 10: (Color online) Blockade voltage VJ, as a function of 
the gate-induced charge in the same SET as in Figs. |H1 and 
The state of the environment is the same as in Figs.l^b) and 
EJb). The dotted (broken) curves denote the potential drop 
across the SET at I = ±0.2 pA (±0.8 pA). 



The relationship between a superconducting SET and 
a single Josephson junction in the low-impedance en- 
vironment has been discussed^ and by examining the 
Hamiltonian, it has been shown that a superconducting 
SET can be viewed as a single Josephson junction with 
gate-tunable Ej. I have found a similar relationship for 
the high-impedance case, which explains the variation of 
\% in Fig. 1101 In the high-impedance environment, not 
only the island charge Q\ — Q2 [— q in Eq. but also 
the total charge (C2Q1 + CxQi)/{Ci + C2) seen from the 
outside contributes to the charging energy^ where Qi is 
the charge on the ith junction. This is why in Fig. [5] the 
Coulomb gap is larger in the high-impedance environ- 
ment. Let us compare Fig. QJb) with the Coulomb gap 
for a single junction in the high-impedance environment, 7 
which is ±e/2C. In this environment, a SET can be 
viewed as a single junction with gate-tunable C or Eq, 
where the minimum of C is ~ 0.25(Ci±C2) and the max- 
imum is ~ 0.5(Ci +C2). I do not know the exact expres- 
sion for the effective Ej of the SET in the high-impedance 
environment, however, the magnitude should be similar 
to the low-impedance case^ < hAo/8e 2 R n . From Tabled 
and C\ + C% = 1 fF, which was obtained from the curve 
fit in Fig. EJ the effective Ej/E c of the SET in Fig. [UJ 



is always smaller than 0.04 (<Cl). Thus, from Eq. JSJ, 
our model predicts that oscillates between ~0.08 and 
~0.16 mV in Fig. The order of magnitude is correct 
and the gate dependence is explained. Note that the no- 
tion of gate- tunable C is a key because at Ej/Ec <C 1, 
Vb depends only on C. 



IV. CONCLUSION 

The transport properties of single-electron transistors 
(SETs) have been studied in various electromagnetic 
environments. In half of the samples, SQUID-array 
leads were employed in order to realize a magnetic-field- 
tunable environment in the superconducting state and a 
high-impedance environment in the normal state. I have 
demonstrated that the effective charging energy of SETs 
in the normal state becomes larger in the high-impedance 
environment than in the low- impedance environment. In 
the superconducting state, the current modulation by the 
gate-induced charge changed from e periodic to 2e peri- 
odic in SETs with SQUID-array leads, as the zero-bias 
resistance R' of the leads was increased to be much larger 
than the quantum resistance Rk = h/e 2 rs 26 kO. This 
change in the periodicity suggests that SQUID arrays 
with sufficiently large R' suppress the contribution of 
quasiparticles to the charge transport. When R' 3> Rk, 
a sharp Coulomb blockade appears in the current -voltage 
characteristics, and the blockade voltage Vb varies de- 
pending on the gate-induced charge. The variation of Vb 
is explained within a model that treats a superconduct- 
ing SET as a single Josephson junction with gate-tunable 
junction capacitance in the high-impedance environment. 
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